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In this article, we calculate the friction between two counter-flowing bosonic and fermionic super¬ 
fluids. In the limit where the boson-boson and boson-fermion interactions can be treated within the 
mean-field approximation, we show that the force can be related to the dynamical structure factor 
of the fermionic component. Finally, we provide asymptotic expressions for weakly and strongly 
attractive fermions and show that the damping rate obeys simple scaling laws close to the critical 
velocity. 
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I. INTRODUCTION 


of the system can be written as 


The onset of frictionless flow in quantum fluids is prob¬ 
ably one of the most intriguing macroscopic manifes¬ 
tations of quantum mechanics. Recent experiments on 
Bose-Fermi superfluid mixtures gave a new twist to this 
old question by probing the critical velocity of a super¬ 
fluid counterflow [1] . When two miscible superfluids flow 
through each-other, Castin et al. suggested a generaliza¬ 
tion of the celebrated Landau criterion where superfluid¬ 
ity is destroyed by the shedding of a pair of elementary 
excitations in the two systems [2] . Later-on, this scenario 
was supported by the study of the lifetime of the quasi¬ 
particles mm or by the calculation of the hydrodynamic 
spectrum [S]. 

The results put forward in [5] were based on heuristic 
arguments and were focusing on the velocity threshold 
above which the counter-flow is damped. In this work, 
we provide a full microscopic treatment of the friction in 
a superfluid counterflow and we determine the explicit 
velocity dependence of the damping force above the crit¬ 
ical velocity. Assuming that the boson-boson and boson- 
fermion interactions can be treated within the mean-held 
approximation, we show that the force can be related to 
the dynamic structure factor of the fermionic superhuid. 
Although the general expression of the structure factor 
of an attractive Fermi gas is not known exactly in the 
crossover between the Bardeen-Cooper-Schrieffer (BCS) 
and molecular Bose-Einstein Condensate (BEC) regimes, 
we provide asymptotic expressions in the limits of weak 
and strong interactions where the fermionic component 
behaves respectively as an ideal Fermi gas and an hydro- 
dynamic Bose-Einstein condensate of dimers [B] . We hnd 
that close to the critical velocity, the force obeys simple 
scaling laws, with exponent depending on the value of 
the fermion-fermion scattering length. 


iL = iLb + iFf + iFbf, (1) 

where iFb,f are the Hamiltonians of the bosons and the 
fermions and iFbf describes the coupling between the two 
species. In the mean-field approximation, we have 

Hbf = ghf j (2) 


where is the field operator of spin a fermions, and 
that of the bosons. Assuming that the BEC can also be 
described within the mean-field approximation, then the 
bosonic field operator can be expanded over the Bogoli- 
ubov modes 


<^(r) = 


rib + 




q 




(3) 

where H is the quantization volume, hkb = mVb, rib is the 
density of bosons and the {uq,Vq) are the usual Bogoli- 
ubov coefficients. The exponential prefactor describes 
the Galilean boost associated with the the motion of the 
condensate (see Appendix]^ . In addition, due to the mo¬ 
tion of the BEC, the energies of the Bogoliubov modes 
are shifted by Doppler effect and their spectrum is given 

by 


Eq.b = E''°lphq-Vb, (4) 

where = hcbq^/Y^-^p^J2 is the Bogoliubov spec¬ 
trum (here Cb is the sound velocity and is the healing 
length) [7]. 

Expanding TFbf to first order in fluctuations of the 
bosonic field, we have 


II. GENERAL SETTING 

We consider a mixture between a fermionic superfluid 
and a weakly interacting Bose-Einstein condensate mov¬ 
ing at velocities V{ and Vt, respectively. The Hamiltonian 


Hm = Hi°'> + + ... (5) 

with 

Em = gurib J ( 6 ) 
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and 


H, 


( 1 ) 

bf 



d^r'^iuq 

q,(T 


Vq)e 


(7) 

being proportional to the number of fermions, it is 
a constant of motion, and can therefore be gauged away. 

and respectively create and destroy Bogoli- 

ubov excitations in the BEC. In the weak coupling limit 
we can expand the damping rate in power of ^bf- To first 
order, only plays a role and using Fermi’s Golden 
Rule the rate of creation of excitations of momentum hq 


in the BEC and in an eigenstate jo;/) of the fermionic 
Hamiltonian H{ is given by 


r(a,q) 


= ^{Uq-VqfS{E^J+Eq,f,)X 

2 

(a/1 J , 


( 8 ) 


where E^j is the energy of the state \af) relatively to 
that of the ground state |0/) of the moving fermionic 
superfluid. Since the creation of a Bogoliubov excitation 
imparts a momentum hq to the BEC, the force acting on 
the bosons can be written as 

F = '^hqV{a,q) (9) 

a,q 

= -Uq)2S'(-q,(10) 

q 


where we introduced the dynamical structure factor of 
the Fermi gas 


S{q,E) = Y,SiE-E^j) 


(a/1 J 


( 11 ) 

Note that the expression for the force is galilean invariant 
and depends only on the relative velocity Vh — V{. Indeed, 
one can show that the dynamical structure factor S' of a 
system moving at velocity Vf is the related to the one of 
a system at rest (denoted as Sgt) by 


Siq,E) = Sstiq,E-hq-Vf) 


( 12 ) 


(this result is a direct consequence of Eq. ( |A12 1 giving 
the energy in the moving frame, and the fact that the 
operator f ^l' 0 cr(r)e*'^ ’' gives a kick hq to the 

system). Inserting this expression in Eq. (10) and using 


the expression for the Bogoliubov spectrum of a moving 
BEC, Eq. 0 we have for the force 


jr = 27 r^ q(uq - VqfS,t(-g, -E^ql - hq ■ F), 


n 


( 13 ) 


with V = Vf — Vb- Finally, using the explicit form of the 
Bogoliubov coefficients, we have 


{Uq - Vqf = 


E 


(0) 

q.b 


(14) 


with eq,b = h‘^q‘^/2mb, hence 

^q.b 




n 


^(O) 
q ^q,b 


5st(-g,-i?S-%-V). (15) 


Eq. (15) is the main result of this paper and shows 


that the friction between the two superfluids is directly 
related to the dynamic structure factor of the fermionic 
component of the system. The dynamical structure fac¬ 
tor describes the response of a given system to a sinu¬ 
soidal perturbation and can be related to two-body cor¬ 
relations. In strongly correlated gases, it was measured 
using Bragg spectroscopy [5] and was used to determine 
the value of Tan’s contact parameter describing the tail 
of the momentum distribution [^. 

The properties of strongly correlated Fermi systems 
are notoriously hard to calculate, and approximate ex¬ 
pression for the dynamic structure factor have been ob¬ 
tained using the approximate methods [lOllllj . However, 
its general properties can be summarized on the sketch of 
Fig.[T} At low momenta, the excitation spectrum is dom¬ 
inated by a phonon branch E = hcfk, with Cf the sound 
velocity of the fermionic superfluid. Since the ground 
state defines the energy zero, Sst vanishes for E < 0. 
However, the presence of a pairing gap implies of a re¬ 
gion where Sst vanishes at positive energies. Landau’s 
critical velocity ul is defined as the highest value of v 
such that Sst{E,q) vanishes for all E < vq (see Fig. [^. 
In the BEC regime of the BEC-BCS crossover, we have 
2 riL = Cf, while in the BCS regime, Landau’s velocity is 
driven by pair breaking excitations and we have ul Cf. 

These general features implies that if for all q we have 
-hq-V < guL, the force is zero. To get some 

damping, there must be values of q for which —E'^I — 
hq ■ V > hqvi^. Using the general Landau argument, we 
get that I V| > ul -I- ra\nq{E^^l/hq) = ul -f Cb. Note that 
this bound sets a necessary, but not sufficient, condition 
on the relative velocity. As a consequence ul + Cb is only 
a lower bound for the real critical velocity. 

Another consequence of this inequality is that the val¬ 
ues of q contributing to the integral are bounded. Indeed, 
from the inequality E < hqv^, we get that 


gCb < \l2x{x F 2), 

where xct, = U — Cb — wl- In other words, when the 
relative velocity gets close to Cb -I- ul, the damping is 
dominated by long wave-length excitations. 

In what follows, we compute the force in the deep BEC 
and BCS limits where analytical expressions can be ob¬ 
tained. 
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FIG. 1: Sketch of the structure factor of an interacting Fermi 
gas: Sat vanishes outside the gray area and A is the excitation 
gap. Dotted line: low momentum phonon branch E = hctq. 
Dashed line: Landau’s critical velocity E = hvi^q below which 
the structure factor vanishes. Note that this sketch corre¬ 
sponds to a weakly interacting system where Landau’s veloc¬ 
ity is hxed by pair breaking, and ul < Cf. In the molecular 
regime of the crossover, vl = ct [21IT2]. 


F 


{271^1%'^ 


d^q 


n 




h,q 


+hq-V)Q{-E^^l-hq-V). 


(16) 

where 0 is Heaviside’s step function. The integral can be 
computed analytically and taking V = H/cb, the force 
becomes 


with 

i(y') 




1 

105v^H 


105H'^ sinh 


1 (6H'® - 39H 


-1 

/4 


- 80 H' 2 -f 8 ) . 


(17) 


(18) 


Close to y' = 1, this expression can be expanded as 

+ (19) 

Even though this scaling is strictly-speaking obtained 
for a non-interacting Fermi gas, it can be extended to the 
case of a BEC moving inside an interacting Fermi gas in 
its normal phase when it can be described within the 
Fermi Liquid Theory. Indeed, close the critical velocity, 
the damping mechanism described involves only low en¬ 
ergy excitations of the system which behave as fermionic 
quasiparticles with renormalized physical parameters. In 
particular, the scaling S{q, E) (x E/q is still valid [Tl] . 


III. BCS LIMIT 


IV. HYDRODYNAMIC LIMIT 


Let’s consider first the deep BCS limit where the 
fermions can be described as an ideal gas. Since there 
is a vanishingly small gap in the fermionic spectrum the 
lower bound derived above is actually saturated and we 
have Vc = Cb. Working close to Cb, the relevant values 
of q much smaller than kp, and we can therefore replace 
the expression of the dynamical structure factor by its 
low-energy approximation (see for instance |13| . and Ap¬ 
pendix 


Sst{q,E) 


rrifflE 

2TT'^hdq' 


this replacement being valid as long as the condition 


kp^ <C 



is fulfilled. 

Replacing the discrete sum over q by an integral, one 
gets 


In the BEC limit of the cross-over, the excitations of 
the fermionic superfluid are dominated by phonons of 
energy Ef,q. In this case, the dynamical structure factor 
can be approximated by a Dirac function and we have in 
this case 


Sat[q,E)=A5{,E-Efl) (20) 

where the normalization constant A is fixed by the f-sum 
rule [Zllli] and is given by 


A = Nf 


'9./ 


E 


( 0 ) ■ 
f,q 


( 21 ) 


In this limit, we obtain for the force 


F = 2:rgl,n,n, ,-|0L«(£») + e'^I j + hg ■ V). 


( 22 ) 

Using the usual Landau argument, we see that the force 
vanishes if the velocity is smaller than a critical velocity 
Vc given by 


Vc = min 
q 


E 


(0) 

q,b 


■E 


(0) ' 
-qj 


hq 


(23) 
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In case of linear dispersion relation, the critical velocity 
is simply 


— Cb + Cf, (24) 

where Cf^b are the sound velocities of the fermions and 
the bosons. 

To calculate explicitly the force, we assume that both 
the bosons and the fermions follow the Bogoliubov dis¬ 
persion relation, an assumption valid in the far-BEC 
limit of the BEC-BCS crossover. Taking ^b.f the healing 
lengths of the two superfluids, we have after performing 
the angular integral 


F 


glfU^nih 


n STTrUfTTlbCbCfE 

riM q^dq ^^\/l + + y -^1 + 




(25) 


with qM given by the condition 

Cb^l + gM^b/2 + Cf^J\F qhi'f 12 = V (26) 
For E ~ I4j qM is vanishingly small, and is given by 


qM — 


V-Vc 

Cb4" + Cf0’ 


hence 


F _ 4g^inhnfh / V - Vg 

n 57rm{mbCbCfVc + Cf0 


5/2 


(27) 


V. CONCLUSION 

In this work we have shown that the damping between 
two bosonic and fermionic superfluids could be related 
to the dynamical structure factor of the fermions. In the 
weakly and strongly attractive limits, we find a power law 
dependence of the damping force vs velocity close to the 
critical velocity. Like previous findings [3], we find that 
the exponent varies across the BEC-BCS crossover and is 
different from the scaling F oc (E' —I)^ found for the drag 
of a particle in a Bose-Einstein condensate [T3]. These 
different scalings highlight the many-body nature of the 
two counter-flowing systems and the role of the collective 
excitations in the damping mechanism. Our prediction 
can be directly tested by measuring the damping of the 
relative oscillations of two superfluids, as presented in [T] 
for a unitary fermions. In particular, comparison with 
Eq. [T^ provides a quantitative test of theoretical pre¬ 
dictions for the structure factor of a strongly correlated 
fermionic superfluid. 
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Appendix A: Bogoliubov Transformation in a 
Moving Frame 


Consider a ensemble of JV particles described by a state 
I'll;). The action of a velocity V boost is described by the 
operator 

fc • , (Al) 

with k = mV/h. In second quantization, this operator 
is represented by 


N 


U = exp / ^ 


i=i 


U = exp 


d^r' ip{r')^ ^(r')k 



(A2) 


Dehne ip'{r) the held operator in the moving frame. 
By dehnition, we must have for 2 ) = t^|'0j=i.2)j 

(i/>(|^'(r)|V' 2 ) = (V'i|^(r’)IV’ 2 ), hence 

(V'lluV {r)U\tp2) = {ipi\0{r)\tlj2) ■ (A3) 

This equation being valid for any \ipi), we must have 
Wif'{r)U = (p{r) hence 

if'(r) = U(p{r)U‘'. (A4) 

(p' can be calculated using the Campbell-Hausdorff for¬ 
mula. Taking U = exp(T), we have indeed 

^'(r) = 0{r) + [f, ip{r)] + ^ [f, [f, v?(r)]] -f ... (A5) 

The commutator [T, (p(r)] is readily calculated and yields 
[T, (^(r)] = —ik ■ ripir). After resumming the Campbell- 
Hausdorff formula, we obtain 

p'ir) = (A6) 


In the moving frame, the BEC is stationary and we can 
therefore decompose p'{r) using a Bogoliubov transform 






LqF‘^"^bq - VqC 


-iqr'fp; 


(A7) 


We thus have in the lab frame 


p{r) = e^^"^p'{r) 




I 






(AS) 
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The eigenenergies of the Bogoliubov excitations of the 
moving Bose-Einstein condensate are obtained by apply¬ 
ing the Galilean boost operator U to the eigenstates |/3)b 
of the condensate at rest. Taking |/3')b = G|/3)b, and 
consider the quantum many-body Hamiltonian 


^ r 

Hb=- J ^(r)d^r 

-b J V{r — r')(p{r)^(p{r')(p{r)^cf'rdPr'. 

(AlO) 


By using Eq. A6), we see readily that if |/3)b is a com¬ 
mon eigenstate of iJb and the momentum and particle 
number operators, then |/3')b is an eigenstate of H\^ for 
the eigenvalue 


Ep = Ejs + hqp -V + Nh 


mV^ 


(All) 


where hqp is the momentum of state |/3)b- Compared to 
the energy Eg of the boosted ground-state Bose-Einstein 
condensate, the energy of the excitation \(3')h is therefore 


E'p-E'o = Ep+ hqp ■ V, (A12) 


and is simply Doppler-shifted. 


Appendix B: Dynamic structure factor of the ideal 
Fermi gas. 




K 


In this appendix, we rederive the dynamical structure 
factor for zero-temperature, non-interacting fermions. A 
more general derivation including finite temperature ef¬ 
fects can be found in [T5] . 

The dynamic structure factor is given by the 
matrix elements of the excitation operator Fq — 
Scr /Expanding the field op¬ 
erators over the plane wave basis as i^air) = 

where ft is a quantization volume, 

yields 


Eg — ^ ‘^k+q,a^k,t7- 
k,a 


i 

1 


K 

q 


K, 


K2 


In other words, the operator Fq kicks one particles of 
the system and gives it a momentum q. The states 
\n) excited by Fq from the ground state |0) correspond 
to particle-hole pairs of momenta k (hole) and k + q 
(particle). The energy of such an excitation is E„ = 
£k+q — £k = /2m -|- ti^k ■ q/m and we can write the 

structure factor as 


S{q,E) = 2 [ 

J k<kF ,\k~\~q\>kF 




2m 


h^kzq 

m 


FIG. 2: Construction of the dynamic structure factor in mo¬ 
mentum space. The shaded area represents the Fermi sphere. 
Energy-momentum conservation fixes the initial kz to a value 
Kz (see text). Since the initial state must lie inside the 
Fermi Sea, we must have k < ki = y/k^ — K/. Moreover, 
the hn al momentum m ust be above the Fermi surface, hence 
K > y/k'j^ + {Kz + qY- The thick solid line represents the mo¬ 
menta contributing to the dynamical structure factor. From 
top to bottom. Upper panel, case (1): Kz -I- q larger than kf. 
The final momentum is always above the Fermi surface. Mid¬ 
dle panel, case 2: Kz -bq is lower than kp, with K 2 < ni. The 
permitted values of k are located in the ring K 2 < k < ki. 
Lower panel, case 3. For K 2 > ki, the contribution to the 
structure factor vanishes. 
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where the z axis is chosen along q and where the factor 
2 conies from the two spin states. 

We perform first the integral over k^- Take K^{q, E) = 

{E — ^^)m/h'^q. The delta-function selects = 
Kz{q, E) and taking k = {kx, ky) we have 


S{q,E) 


2m f (pKfl 

^^<7 J\K^\<kF,\k+q\>kF (27^)^ 


(Bl) 


The condition \Kz\ < kp selects a band in the {q,E) 
plane. Indeed, we must have 


h^kpq ^ ^ fi^kpq 

m ~ \ 2m J ~ m 

The constraint kz = Kz selects a plane in initial- 
momentum space. Taking into account the fact that the 
initial momenta must lie within the Fermi Sea, we see 
that the transverse momentum k lies within a disk of ra¬ 
dius Ki = \Jk‘p — Kl- The set of allowed values for k is 
further constrained by Pauli Blocking on the final states. 
Indeed, the final momentum must lie outside the Fermi 
Sea and as such it implies the condition |fc -f q| > kp- 
We must then consider three cases (see also Fig. : 

1 . If Kz -I- <7 is larger than kp , then for any k such that 
kz = Kz, the final momentum fe -|- g is outside the 
Fermi Sea. The constraint on the final state is thus 
always satisfied and does not have any effect on the 
calculation of the integral over k, which yields the 
area of the disk of radius . We thus have 


For E < 0, the dynamic structure factor vanishes, 
as expected since the Fermi Sea is the ground state 
of the system. 



FIG. 3: Dynamical structure factor of the ideal gas. The 
shaded area corresponds to the region where the structure 
factor takes hnite values. The solid lines corresponds to con¬ 
dition — kp < Kz < kp. The dashed marks the limit be¬ 
tween cases (1) and (2) and is associated with the condition 
Kz+q = kp. 


S{q,E) 


mflui 


Finally, let’s consider some simple limiting cases. For 
q > 2kp, we are in case (1) (see Fig. [^, and S{q,E) is 
(B2) therefore given by 


2. For Kz + q < kp, final momenta corresponding 
to small values of k (more precisely smaller than 
K 2 = y/kp — {Kz + qY) lie inside the Fermi Sea 
and are therefore forbidden by Pauli Principle. We 
have then 


S{q,E) 


mfl{Kf — K2) 


m^QE 

2p^h^q 


(B3) 


Take kp = rif, where ni is the total density of 
fermions, this expression can be recast as 


S{q,E) 


mfl 2 (p k?q^\^ 

Ep'^h'^q ^ ti^q^ \ 2m ) 


(B5) 


In this regime, S is peaked around the energy E = Eg = 
h?q^/2m with a half-width IS.E = h?kpq/m. We note 
that /S.E/£q = 2kp/q ^ 0 for q ^ 00 . In other words, we 
recover at large momenta the dynamic structure factor 
of a free particle. 

For low values of q, the range of validity of case (1) 
is vanishingly small and the dynamic structure factor is 
mostly given by case (2). We therefore have 


S{q,E) 


3N{ /kp\ E 


(B4) 


S{q,E) 


3N{ /kp\ E 


(B6) 


with TVf = rif^l the total number of fermions. 

3. Eq. |B3| is valid as long as Ki > ^ 2 - Graphically, 
case (2) ends when Kz + q = —Kz, ie for E = 0. 


as long as 0 < E < 2Epq/kp, and ^ = 0 otherwise. 

One readily checks that both approximate expressions 
(B6) and (B5) satisfy the f-sum rule. 
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